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Abstract 

This paper investigates the evolution of the dust-to-metal ratio in galaxies based on a simple evolution 
model for the amount of metal and dust with infall. We take into account grain formation in stellar mass- 
loss gas, grain growth by the accretion of metallic atoms in a cold dense cloud, and grain destruction by SNc 
shocks. Especially, we propose that the accretion efficiency is independent of the star-formation history. 
This predicts various evolutionary tracks in the metallicity (Z)-dust-to-gas ratio (T>) plane depending on 
the star-formation history. In this framework, the observed linear Z—T> relation of nearby spiral galaxies 
can be interpreted as a sequence of a constant galactic age. We emphasize that an observational study 
of the Z—T> relation of galaxies at z ~ 1 is very useful to constrain the efficiencies of dust growth and 
destruction. We also suggest that the Lyman break galaxies at z ~ 3 have a very low dust-to-metal ratio, 
typically ^0.1. Although the effect of infall on the evolutionary tracks in the Z-T> plane is quite small, 
the dispersion of the infall rate can disturb the Z-T> relation with a constant galactic age. 
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1. Introduction 

Because dust exists everywhere, stellar light is dimmed 
by the interstellar dust, strong far-infrared radiation from 
dust is observed, and the interstellar gas is heated/cooled 
by some processes involving dust. Moreover, comets and 
planets are made from dust, and various molecules are 
formed on dust. Therefore, dust is one of the most im- 
portant constituents of galaxies and the universe (e.g., 
Whittct 2003). 

In the context of research on the galaxy evolution, 
knowledge of the dust amount in galaxies is essential to 
extract the physical properties of galaxies from observa- 
tions. This is because radiation from stars is always ex- 
tinguished to some extent by internal dust attenuation. 
Since it is difficult to estimate the precise amount of at- 
tenuation from the observed photometric data (e.g., Buat, 
et al. 2002), any other information about the amount of 
dust is very useful. One may think that the metallicity is 
a good tracer of the amount of dust. This is true if the 
dust-to-metal ratio can be regarded as universal. 

Dust grains arc made from metal elements, which 
are injected into the interstellar medium (ISM) by stel- 
lar mass-loss at the end of, or during, stellar evolu- 
tion. The evolution of the amount of metal in a galaxy 
is traced as a consequence of the star formation ac- 
tivity of the galaxy. To date, such metallicity evolu- 
tion has been studied very well since classical work in 
the 1970s (e.g., Audouze, Tinsley 1976; Tinsley 1980 
for reviews). The evolution of the amount of dust can 
be traced by an extension of the evolution of metal- 
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licity if we take account of the formation, destruction, 
and growth processes of dust grains in the ISM. After 
the pioneering work by Dwek, Scalo (1980), some au- 
thors have examined this subject (Wang 1991; Liscnfcld, 
Ferrara 1998; Dwek 1998; Hirashita 1999a,1999b,2000a; 
Edmunds 2001;Hirashita, et al. 2002a,2002b;Hirashita, 
Ferrara 2002; Morgan, Edmunds 2003). 

For several nearby spiral galaxies, there is a global linear 
correlation between the metallicity (Z) and the dust-to- 
gas ratio (T>) suggested by Issa, et al. (1990) observation- 
ally. This means that the dust-to-metal ratio is nearly 
constant for these spiral galaxies. A simple model de- 
veloped by Hirashita (1999a) reproduces this global Z-T> 
relation well. Edmunds (2001) argues that the dust-to- 
metal ratio does not change significantly during the entire 
galactic evolution based on his simple model. In short, the 
global Z~T> relation can be interpreted as an evolutionary 
sequence by the model of Edmunds (2001): galaxies evolve 
from a point of lower Z and T> to another point of higher 
Z and T> with keeping a constant dust-to-metal ratio. 

On the other hand, the Z-T> relation for dwarf galaxies 
is not so tight; the dispersion reaches about an order of 
magnitude. Liscnfeld, Ferrara (1998) propose that out- 
flows triggered by star-formation activity can cause such 
a large dispersion, whereas Hirashita et al. (2002b) show 
that the dispersion can result from a large variety of the ef- 
ficiency of dust destruction by supernovae (SNe) shocks, 
even if there are no outflows; in their scenario the va- 
riety of the efficiency is caused by an intermittent star- 
formation history. In any case, the Z-T> relation (or dust- 
to-metal ratio) is not universal, at least for dwarf galaxies. 

What determines the dust-to-metal ratio? It is a very 
important question. Metal elements in the high-z universe 
have already been detected (e.g., Pettini, et al. 1997). To 
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know the amount of dust (or extinction) , we need to know 
what fraction of metals is condensed into dust. Recently, 
there are several models which can be used to calculate 
the spectral energy distribution of galaxies by solving the 
radiative transfer in the dusty ISM coupled with the chem- 
ical evolution model (e.g., Silva, et al. 1998; Takagi, et 
al. 2003). In these models, a dust-to-metal ratio, such as 
that in the Galaxy, is assumed. Is this really correct for 
dwarf galaxies, starburst galaxies, and high-z galaxies? 

In this paper, we investigate the evolution of the dust- 
to-metal ratio in galaxies by using a simple model for the 
evolution of the amount of dust and metal. Although 
Dwek (1998) presents the most detailed evolution model 
for the amount of metal and dust abundance for our 
Galaxy, we need not to trace so the detailed abundance of 
metal and the composition of dust. The aim of this paper 
is to discuss on a global (or mean) dust-to-metal ratio in 
a galaxy. Therefore, a simple treatment, as adopted in 
Hirashita's and Edmunds's work, is sufficient. 

In the next section, we describe the model used in this 
paper. In the model, we take account of gas (and also 
metal) infall in order to discuss the time evolution on a real 
time-scale. Moreover, we reconsider the efficiency of the 
grain-growth process in the ISM. We then show in section 
3 that the observed constancy of the dust-to-metal ratio 
in nearby spiral galaxies can be interpreted as a sequence 
of a similar galactic age. The achieved conclusions in this 
paper are summarized in the last section. 

2. Model 

Here, we describe the model used in this paper. First 
a set of equations and meanings of notations are summa- 
rized. We next explain the considered processes of dust 
formation, growth, and destruction in the ISM. Then, a 
set of parameters suitable for our Galaxy is estimated. 

2.1. Equations 

First of all, the basic assumptions adopted here are sum- 
marized as follows: (1) the ISM is treated as one zone, but 
two components. One is a cold dense cloud, and the other 
is a warm diffuse medium. The mass fraction of the cold 
cloud is fixed as a constant throughout the whole evolu- 
tion of a galaxy. (2) Stars are formed only in the cold 
dense cloud following the Schmidt law of index 1. (3) The 
instantaneous recycling approximation is adopted. (4) 
Produced metal and dust are quickly mixed well in the 
ISM. (5) Dust consists of two components. One is the re- 
fractory core, and the other is the volatile mantle, which 
can exist only in a cold cloud. In other words, the mantle 
dust is destroyed quickly if it goes out of the cloud. (6) 
No outflow of the galactic scale. 

Under the above assumptions, the time evolutions of the 
gas mass (M gas ), the metal mass (Mz), the core dust mass 
(Mcor), and the mantle dust mass (M man ) are governed by 
the following equations: 



Table 1. Notations and standard set. 
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Parameter 


Value 


Meaning 


T"SF 


5 Gyr 


star formation time-scale 


An 


0.5 


infall parameter (= rap/im) 


R 


0.299 


returned mass fraction 


yz 


0.025 


stellar yield 


Zin 


0.01Z Q 


metallicity in infall gas 


fc 


0.1 


condensation efficiency 


PSN 


10 


SNe parameter (= tsf/tsn) 


7acc,0 


10 8 yr 


accretion growth time-scale 


V 


0.5 


cold gas fraction 




0.5 


refractory probability 



dt r S F 



T~SF 



Z in F, 
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+ fo[R 



,yMz 

T~SF 



yz- 



T)M ( 



TSF 



Mrr 



+ (1-0 



M cor (l - 6 cl ) 
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(4) 



dt t sf 

where Tgp, TgN, and r acc are the time scales of the star 
formation, dust destruction by SNe, and dust growth by 
metal accretion in the cold cloud, respectively. The term 
yyAfgas/TsF is the star-formation rate, where rj is the mass 
fraction of the cold gas. The notations are summarized in 
table 1 and are described in the rest of this subsection. In 
section 3, we will solve equations (1) — (4) combined with 
the initial condition, M gas = Mz = M COI = M man = at 
t = 0. 

We consider gas and metal infall from a halo, which 
can resolve the so-called G-dwarf problem of the stellar 
metallicity distribution in the Galaxy (e.g., Twarog 1980). 
Here, the gas infall rate is denoted by F, which is assumed 
to be the following analytic function: 

F=^exp[-f/T in ], (5) 

Tin 

where r ln is the time-scale of infall and Mq is the total 
gas mass that can fall into the disk from the halo i.e., 
J Fdt = Mq. The metallicity in the infalling gas is denoted 
as Z[yi • 

For the chemical evolution part [equations (1) and (2)], 
we need the returned mass fraction, R, and the stellar 
yield, yz- The returned fraction is the mass fraction of 
returned gas from stars into the ISM per unit mass of star 
formation, defined as 

/>m up 

R= (m — w m )4>{m)dm , (6) 

J ni t 

where w m is the remnant mass of a star with mass m, m t 
is the present turn-off mass, and <fi(m) is the initial mass 
function (IMF), normalized as 



m<f>(m)dm = 1 , 



(7) 
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where m\ ow and m up are the lower and upper cutoff masses 
of the IMF. When we assume the Salpeter IMF with 
Wlow = 0-IMq, m up = 100M Q , and m t = 1M , and as- 
sume the remnant mass function in Pagel (1997), we ob- 
tain R = 0.299. The stellar yield is the mass fraction of 
the metal elements newly produced and ejected into the 
ISM per unit star formation, and defined as 



Vz 



mpz(m)(j)(m)dm , 



(8) 



where pz{m) is the fractional mass of metal elements 
ejected into the ISM by a star of mass m. Assuming the 
Salpeter IMF with the same characteristic masses above 
and pz{ m ) m Pagel (1997), we obtain yz = 0.025. 

In the part concerning the evolution of the amount of 
dust [equations (3) and (4)], f c is the condensation effi- 
ciency in the stellar mass-loss gas, and S c \ is the dust-to- 
metal ratio in a cold cloud, i.e., 



M„ 



M z ' r)M z ' (9) 

Finally, £ is the refractory probability of metal accretion. 
The refractory core may grow by an accretion process in 
the cold cloud as well as the volatile mantle. This pa- 
rameter is defined as the fraction of accreted metal mass 
which contributes to the growth of the refractory core. 

2.2. Grain Formation, Evolution, and Destruction 

We first consider the condensation efficiency in the stel- 
lar outflow, f c . The efficiency may depend on the stellar 
mass; different efficiencies between that in mass ejection 
by SNe and that in stellar wind from evolved stars are 
possible. In the framework presented here, we do not 
distinguish these two types of stellar mass-loss because 
of the instantaneous recycling approximation. Thus, the 
condensation efficiency given in this paper is an effective 
one, including both types of mass loss and averaged over 
the stellar mass by the assumed IMF. Moreover, the ef- 
ficiency can evolve with the stellar metallicity. Indeed, 
Morgan, Edmunds (2003) suggest such an evolution, but 
the evolution is not very significant, fortunately. This is 
because the metal amounts in stellar outflows are often 
dominated by the newly produced metal. Therefore, we 
neglect the metallicity dependence of / c , and assume a 
constant f c . 

The "stardust" injection rate is f c (RZ + yz)r]M geis /TsF 
[i.e., third term in equation (3)], where Z is the metallicity 
(= Mz/M gas ). On the other hand, the stellar mass- loss 
rate is Rr]M sa3 /TsF- Again, we note that the injection rate 
and the mass-loss rate include all types of stellar outflows. 
Thus, the dust-to-gas ratio in the mass-loss gas (2?*) is 



f c (RZ + y z ) 
R 



, yz 
Jc r ■ 



(10) 



The approximation in the last term is valid when Z <IC 
4Z Q for R = 0.299, y z = 0.025, and the solar metallicity 
Zq = 0.02. Observationally, 23* ~ 0.01 is suggested as 
the present-day Galactic value (Gehrz 1989), although the 
uncertainty is still large. Therefore, we obtain / c ~ 0.1, 



which agrees with values given by Hirashita (1999a, 1999b) 
and Edmunds (2001). 

Next, we discuss the grain growth rate in the cold cloud. 
The dust mass increases if the metallic atoms accrete on 
the pre-existing grain cores. The refractory core, itself, 
may grow by accretion as well as the volatile mantle. If the 
accreted metal elements are refractory, such as Si, Fe, or 
Ti, the refractory core grows. On the other hand, elements 
such as C, N, and O produce a volatile mantle (maybe 
ice molecules) on the refractory core when these elements 
accrete onto the grain surface. Actually, the accretion 
rates of various elements are different from each other. 
To avoid such difficulty, we parameterize the situation by 
the refractory probability, £. This is defined as the mass 
fraction of accreted metal contributing to the growth of 
the refractory core. 

The total (corc+mantlc) dust mass growth rate by 
metal accretion can be expressed as 



dAL 



dt 



= N? or ira 2 szp c z(vz): 



(11) 



where N^ r 1S the number of refractory cores in the cold 
clouds, a is the typical grain radius, sz is the mean stick- 
ing coefficient of metal elements, p c z is the mass volume 
density of gaseous metallic atoms in the cold clouds, and 
{vz) is the mean velocity of metallic atoms. We assumed a 
spherical grain for simplicity. Since N^ ov = r)M CO v/rnc\ and 
p c z = p gas Z(l — ^ci)j where is the typical grain mass 
and p gas is the volume mass density of the cold cloud gas, 
the dust mass growth rate is reduced to 



M cor (l-<5 cl ) 



dt 

where we 



(12) 



define the accretion time scale as r acc = 
m ( i/riTra 2 szPg as {vz)Z. We then obtain the fourth term 
in equation (3) and the second term in equation (4) by 
using the refractory probability, £. The parameter is de- 
termined later so as to reproduce the Galactic metallicity 
and dust-to-gas ratio. 

The accretion time-scale, r acc , is not constant during 
galaxy evolution, but depends on the metallicity. Here, 
we normalize it as 



(13) 



and we estimate the normalization, r acCj o, which means 
the accretion time-scale at the solar metallicity (i.e., 
the current value in the Galaxy). According to 
Hirashita (2000b), who estimated the grain growth time- 
scale by taking into account the mass function of the 
molecular clouds in the Galaxy, r acc ^ ~ 77" 1 5 x 10 7 yr. 
Although the cold gas fraction, r\, is rather uncertain, r\ 
may not deviate significantly from 0.5. Thus, we always 
assume T acCi o ~ 10 8 yr, except when indicated otherwise. 

We next consider the dust-destruction rate. There are 
two processes for destruction; one is incorporation by star 
formation; the other is evaporation and sputtering in SNe 
remnants. Here, we assume that star formation occurs 
only in a cold cloud. Since the dust-to-gas ratio in a cold 
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cloud is {rjMcox + M ma , n )/r)M sas , and the star-formation 
rate is r)M gas /TSF, we obtain the first terms in equations 
(3) and (4) as the destruction rate by star formation of 
core and mantle, respectively. Here, wc have assumed the 
destruction efficiency by star formation to be unity. 

The SNe destroy only the core dust because the mantle 
dust is already destroyed in the SNe remnants, which are 
not cold gas. In fact, the SNe affect the amount of mantle 
dust through the cold gas fraction, which is reduced by 
the SNe. However, we assume a constant fraction of cold 
gas for simplicity. Since the dust-to-gas ratio in a warm 
medium is M cor /M gas , the destruction rate by SNe is 



dM d 
dt 



Mr, 



SN 



M„ 



-eM SNR 7, 



(14) 



where e is the destruction efficiency, Mg^R is the gas mass 
shocked by one SN, and 7 is the occurrence rate of SNe. 
We do not distinguish type la and type II SNe. Thus, 
the SNe-rate 7 means the sum of the occurrence rates of 
both type la and type II. If we define the time-scale of the 
destruction by SNe as r SN = M gas /eM SNR 7 (McKee 1989), 
the second term in equation (3) is obtained. 

Dwek (1998) proposed TsN/T aC c,o — 2 (constant). It was 
also assumed in Hirashita (1999a, b). Here, we reconsider 
this argument. When the star- formation rate is denoted 
as ip(= t]M S8iS /t sf ), we find r SN /r SF = /(^cMsnr) • 

On the other hand, we have seen that r acCi o ocrj^ 1 . Thus, 
we find TsN/T a cc,o (V'/7) t sf- We may consider that the 
ratio of ip and 7 is roughly constant. Therefore, the ra- 
tio rsN/T aC c,o is proportional to tsf- In other words, the 
ratio is not constant, as assumed by Dwek (1998), but 
dependent on the star-formation history (i.e., Tsf)- This 
provides us with a new interpretation of the Z-T> relation 
for nearby spiral galaxies presented in the next section. It 
is worth noting that we do not need the assumption of a 
constant 77 in the above statement. 

Although the value of r\ is required when we esti- 
mate r acc .o and tsf, estimating 77 is a rather complex 
task. After stars are formed and SNe occurs, the cold 
gas is heated by the SNe shocks, so that 77 decreases. 
When 77 decreases, the star- formation rate is suppressed, 
and then the SNe occurrence rate decreases, so that 77 
increases. Thus, the star- formation rate (oc r^p) and 
77 are coupled with each other non-linearly (Ikeuchi, 
Tomita 1983; Hirashita 2000a). To avoid the complex 
task of solving a set of non-linear equations to determine 
77, we assume a constant 77 as an average value during 
the galaxy evolution time-scale for simplicity. For a con- 
stant 77, the ratio of t§n to tsf becomes a constant. Since 
McKee (1989) estimates tsn — 4 x 10 s yr for the Galaxy 
and tsf ~5x 10 9 yr is suitable for the Galaxy, we obtain 
tsf/tsn ~ 10. Combining T acCj0 — 10 s yr, we find 



tsn 

Tacc,0 



tsf 
1 Gyr 



(15) 



In future work, the fraction 77 should be determined sclf- 
consistcntly in the model. 




Fig. 1. Time evolution of the metallicity, dust-to-gas ratio, 
and dust-to-metal ratio. Bottom : the solid curves indicate 
the metallicity and the dashed curves indicate the dust-to-gas 
ratio. The vertical axis is normalized by the present-day 
Galactic value, which is indicated by the dotted straight line. 
The thick curves are the standard case (table 1), the top thin 
(solid and dashed) curves are the case of r\ = 0.8 and £ = 0, 
and the bottom thin curves are the case of r) = 0.2 and f = 1, 
where 77 is the cold gas fraction and £ is the refractory prob- 
ability in the accretion process. Top : the solid, dashed, and 
dash-dotted curves are the standard case, r\ = 0.8 and 5 = 
case, and 77 = 0.2 and £ = 1 case, respectively. The dotted 
straight line indicates the present-day Galactic value. 



2.3. Standard Case (Milky Way) 

Here, wc estimate a set of parameters suitable for our 
Galaxy. First, the time-scales of the star formation and 
infall are discussed. Takeuchi, Hirashita (2000) have ex- 
amined the star-formation history of the Galaxy by us- 
ing a chemical evolution model with infall similar to that 
used in this paper. According to them, a set of these 
time-scales between (rsF,Tin) = (6,23) and (11,12) in Gyr 
unit can reproduce the observed star-formation history of 
the galactic disk, the age-metallicity relation of stars in 
the solar neighborhood, and the metallicity distribution 
of G-dwarfs. Noting that their star-formation rate and 
that in this paper are different in the treatment of the 
cold-gas fraction (77), we adopt t$f = 5 Gyr and r- m = 10 
Gyr. Combining these time-scales with R = 0.299 and 
yz = 0.025 (sec subsection 2.1), and assuming a very low 
metallicity in the infall gas (for example Z- m = O.OIZq), 
we find that the present-day (age ~ 12 Gyr) metallicity 
becomes around Solar metallicity (see figure 1 bottom 
panel), which ensures the validity of these time-scales. 

For the convenience of later discussions, we define the 
infall parameter, /3- in , as rsF/Ti n - Thus, we adopt f3- la = 
0.5 for our Galaxy. Since an infall parameter suitable 
for other galaxies with a different star-formation history 
(i.e., tsf) is rather unknown, we assume (3- m = 0.5 for all 
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galaxies. The effect of the other choices of f3- nl is discussed 
in subsection 3.2. 

We use values estimated in subsection 2.2 for parame- 
ters f c , /3 SN ee t sf /t sn . and r acCj0 . Strictly speaking, r acCj0 
depends on 77, as shown in subsection 2.2. Nevertheless, 
we adopt a constant r acc = 10 8 yr for simplicity because 
the uncertainty of r acCi o is large and 77 may be roughly 
-0.5. 

We now modulate two parameters of 77 (cold gas frac- 
tion) and £ (refractory probability in accretion) so as 
to reproduce the present-day dust-to-mctal ratio in the 
Galaxy. Adopting the solar metallicity (0.02) and the 
dust-to-gas ratio in the solar neighborhood (6 x 10 -3 ; 
Spitzer 1978) as the present-day Z and V in the Galaxy, 
we find the present-day dust-to-metal ratio, 8 = (M cor + 
-Mman)/Afgas = 0.3. In figure 1, we show the time evolution 
of the metallicity, the dust-to-gas ratio, and the dust-to- 
metal ratio for three sets of 77 and £: (1) 77 = £ = 0.5, 
(2) 77 = 0.8 and £ = 0, and (3) r) = 0.2 and £ = 1. All of 
these cases give the calculated 8 ~ 0.3 at the galactic age 
of 12-15 Gyr (top panel in figure 1). From the bottom 
panel of figure 1, we find that the case of 77 = £ = 0.5 is the 
best. Thus, we call this case the 'standard case', whose 
parameters are summarized in table 1. 

Although the case of 77 = 0.2 and £ = 1 is not suitable 
for the Galaxy, because of the too small metallicity, the 
case of 77 = 0.8 and £ = is not so bad. Since the time 
evolution of the dust-to-metal ratio in this case is very 
similar to that of the standard case, our results obtained 
later are not altered if we choose this case as the standard 
case. We note here that the parameter set given in table 
1, as the standard set, is not a unique solution for the 
Galaxy. There may be other sets which can reproduce all 
observational constraints for the Galaxy. However, we do 
not have enough data for determining a unique solution, 
even for the Galaxy. For example, we can determine £ 
if we know the mass ratio of the grain core and mantle, 
which is still unknown. In future work, we may determine 
these parameters more precisely. 

It is important that the dust-to-metal ratio in young 
galaxies can be much smaller than that in the present-day 
Galaxy, as shown in the top panel of figure 1. Therefore, 
a galaxy does not evolve with a constant dust-to-metal 
ratio. That is, we cannot consider the Z—V relation of 
Issa et al. (1990) to be an evolutionary sequence, like 
Edmunds (2001). We discuss what the Z-V relation im- 
plies in the next section. 

3. Results and Discussions 

3.1. Z-V Relation 

Here, we compare the observed Z-V relation with our 
theoretical ones. First, the observational data used are 
explained. Although a large and systematic data set of Z 
and V of various types of galaxies is required, we cannot 
find such a data set. Thus, we use the data presented by 
Issa ct al. (1990), whose sample galaxies are five nearby 
spiral galaxies (the Galaxy, M 31, M 33, M 51, and M 101) 
and both Magellanic clouds. Lisenfeld, Ferrara (1998) 
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Fig. 2. Metallicity (Z)-dust-to-gas ratio (T>) relation. The 
vertical and horizontal axes are normalized by Zq and T> of 
the Galaxy, respectively. The open circles are observed data of 
Z and T> for five nearby spiral galaxies and Magellanic clouds 
taken from Issa et al. (1990), except for M 31 whose T> is taken 
from Inoue (2001). The filled circle is M 31 taken from Issa 
et al. (1990). The solid curves are theoretical Z—T> relations; 
the adopted star formation time-scale tsf = 1, 2, 5, 10, 20, 
50, and 100 Gyr from top to bottom. Other parameters are 
the same as those in table 1. The dash-dotted curves indicate 
the iso-galactic age; 1, 5, 10, and 15 Gyr from left to right. 
That is, the thin, medium, and thick shaded areas indicate 
that the galactic age is in 1—5 Gyr, 5-10 Gyr, and 10-15 Gyr, 
respectively. The dashed straight line represents the sequence 
of the same dust-to-metal ratio as the Galaxy. 

and Hirashita et al. (2002b) compiled such data of dwarf 
galaxies. However, since their Vs of dwarf galaxies are 
based only on IRAS data, there may be systematic un- 
derestimations due to the cold dust not detected by IRAS 
(e.g., Popescu, et al. 2002). Thus, we decided to use only 
Issa's data. For M 31, the dust-to-gas ratio estimated in 
Inoue (2001) from the data obtained by ISO was also used. 
It is very useful to examine the Z-V relation for a large 
sample of galaxies in future work. 

In figure 2, we display a comparison between the ob- 
served data points and the theoretical predictions with 
various star-formation histories (i.e., rgp). The adopted 
tsfs are 1, 2, 5, 10, 20, 50, and 100 Gyr from top to 
bottom. They are suitable for early- to late-type galaxies 
from top to bottom. In our model, a different evolution- 
ary path in the Z-V plane is predicted depending on tsf- 
This is because TSNAacc,o is not constant, but dependent 
on tsf, as discussed in subsection 2.2 [equation (15)]. In 
other words, the accretion time-scale does not depend on 
TsF) whereas the SNe destruction time-scale depends on it. 
The separation of these evolutionary sequences is caused 
by the different efficiency of the dust destruction by SNe 
relative to the accretion growth of dust. A more rapid 
star formation causes a higher occurrence rate of SNe, 
and results in a more efficient destruction of dust. 

How can we understand the observed constancy of the 
dust-to-mctal ratio in nearby spiral galaxies within the 
framework of our model? The answer is that the sequence 
of the dust-to-metal ratio is the sequence of a constant 
galactic age. As shown in figure 2, these spiral galaxies 



G 



A.K.Inouc 




0.01 0.1 



dust-to-gas ratio 

Fig. 3. Predicted redshift evolution of Z— T> relation. The 
solid curves are the predicted Z—T> relation at z = 3, 2, 1, 0.5, 
0.1, and (local) from left to right, as indicated in the panel. 
The dotted straight line represents the sequence of the same 
dust-to-metal ratio as the Galaxy. The vertical and horizontal 
axes are the same mean as those of figure 2. The present-day 
galactic age is assumed to be 12 Gyr. The adopted cosmology 
is Q M = 0.3, tt A = 0.7, and H = 70 km s" 1 Mpc" 1 . The 
vertical and horizontal axes are normalized by Zq and T> of 
the Galaxy, respectively. 

are plotted in (or around) the area with an age of 10-15 
Gyr (thick shaded area). Therefore, we conclude that the 
observed linear Z-T> relation for the nearby spiral galaxies 
is not an evolutionary sequence with a constant dust-to- 
metal ratio, but a sequence of a similar galactic age with 
various star- formation histories (i.e., from early- to late- 
type). This is the main conclusion of this paper. 

As shown in figure 1, our model predicts a time evo- 
lution of the dust-to-metal ratio: a lower dust-to-metal 
ratio in younger galaxies. Thus, the Z-T> relation is ex- 
pected to evolve with the galactic age or redshift. We 
show such redshift evolution of the Z-T> relation in figure 
3, where we assume the galactic age at the present-day 
to be 12 Gyr. We predict a deviation from the local Z- 
V relation at z > 0.5. Therefore, the investigation of the 
Z~T> relation in a high-z will be a good test against our 
model. Observations with, for example, ASTRO-F and 
SIRTF will be very useful. 

Finally, we comment on the Magellanic clouds. They 
are plotted in a somewhat younger area in figure 2 than 
other spiral galaxies. They are found in the area of an age 
of 5-10 Gyr (medium shaded area). The definition of the 
galactic age in the model is the time from the start of the 
main gas infall and star formation because our treatment 
is one-zone. If a minor star formation exists before the 
main one, we underestimate the galactic age. Thus, the 
younger ages of Magellanic clouds may be due to such an 
effect. 

3.2. Effect of SNe, Accretion, and Infall Time-Scales 

In this subsection we discuss the effect on our inter- 
pretation of the local Z-T> relation when we change tsn , 
Tacc.Oj an d T; n . Edmunds (2001) argues that the accre- 
tion time scale is much shorter than that adopted here, 
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Fig. 4. Same as figure 2, but T aC c,o = 10 7 yr and rj = 0.3 are 
adopted. The thin, medium, and thick shaded areas indicate 
the same areas as in figure 2 for a comparison. The verti- 
cal and horizontal axes are normalized by Zq and T> of the 
Galaxy, respectively. 

and that the dust destruction by SNe is not efficient. 
Indeed, he assumes the sticking probability on the grain 
surface in the accretion process of the metallic atoms to 
be unity, whereas our estimate is based on the value of 
Hirashita (2000b), who assumes a sticking probability of 
0.1. The accretion time scale becomes an order of magni- 
tude shorter than that in table 1 if we assume a sticking 
probability of unity. 

In figure 4, the high accretion cases are displayed. All 
curves and points, except for the shaded areas, have the 
same meanings as in figure 2, but r acc .o = 10 7 yr and r\ = 0.3 
are adopted. The value of r\ is determined so as to repro- 
duce the parameters of the Galaxy, as done in subsection 
2.3. The shaded areas cover the same areas as figure 2 
for a comparison. We find that the variation of the Z— 

V relations becomes small and the evolutionary tracks of 
various tsf are almost superposed on the dashed straight 
line, which is a sequence of the same dust-to-metal ratio 
as the Galaxy if the galactic age is larger than about 5 
Gyr. This means that the local Z-D relation of spiral 
galaxies can be interpreted as the evolutionary sequence. 
Of course, we can interpret it as a sequence of similar 
galactic age, although sequences of age larger than about 
5 Gyr cannot be distinguished from each other. 

In figure 5, we display the case of no SNe destruction. 
All curves and points, except for the shaded areas, have 
the same meanings as in figure 2, but /3sn = 0, rj = 0.3, and 
£ = are adopted. The values of 77 and £ are adjusted so as 
to reproduce the Galactic parameters. The shaded areas 
indicate the same areas as in figure 2 for a comparison. We 
find a similar trend to the high accretion case, a tighter Z— 

V relation. We can again interpret the local Z-D relation 
as the evolutionary sequence. 

Although the data point of the Large Magellanic Cloud 
cannot be covered by the no SNe destruction and high 
accretion unless tsf < 1 Gyr, this may not be so critical 
because the uncertainty of the data point is likely to be 
larger than its error-bars shown in figures 2, 4, and 5. 
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Fig. 5. Same as figure 2, but /3gN = Tgp/rgN = 0, r) = 0.3, 
and 5 = are adopted. The thin, medium, and thick shaded 
areas indicate the same areas as in figure 2 for a comparison. 
The vertical and horizontal axes are normalized by Zq and 
T> of the Galaxy, respectively. 
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Fig. 6. Time evolution of the dust-to-metal ratio. The ver- 
tical axis is normalized by the present-day value. The solid, 
dotted, and dashed curves are the standard (table 1), no SNc 
destruction, and high accretion cases, respectively. For all 
curves, tsf = 5 Gyr and the present-day galactic age of 12 Gyr 
are assumed. The adopted cosmology is Qjy[ = 0.3, Ca = 0.7, 
and H = 70 km s _1 Mpc" 1 . 

Therefore, we cannot distinguish the standard, no SNc 
destruction, and the high-accretion cases only by the local 
Z-D relation. However, fortunately, these three cases can 
be distinguished if we have information about the Z—T> 
relation of younger galaxies. 

In figure 6, we show the dust-to-metal ratios of the stan- 
dard, no SNe destruction, and high accretion as a function 
of the lookback time and redshift. A present-day galactic 
age of 12 Gyr and tsf = 5 Gyr are assumed. From this fig- 
ure, we find that the dust-to-metal ratios of the standard 
and no SNe destruction cases at z = 1 become 30% and 
60% of the local value, respectively, whereas that of the 
high-accretion case does not change significantly. Thus, 
we can obtain further constraints on t§n and T acC: o by us- 
ing observations of the Z-T> relation at z ~ 1 (or local 
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Fig. 7. Effect of infall on the Z— T> relation. The solid and 
dotted curves are the /3; n = 0.5 and no infall cases, respectively. 
The adopted rgp = 1, 2, 5, and 10 Gyr from top to bottom. 
The vertical and horizontal axes are normalized by Zq and 
T> of the Galaxy, respectively. 

galaxies with an age of several Gyr). 

The effect of infall on the Z-V relation is discussed. In 
Figure 7, we show the comparison between the standard 
cases and the no infall cases. Clearly we find the effect of 
infall is negligible in the plane of Z and V. However, the 
infall parameter j3- ln affects the speed of the chemical en- 
richment. Thus, a galaxy without infall rapidly evolves in 
the Z-T> plane along the track determined from the com- 
bination of parameters, r acCj o, /3sn, ^7, and £. Therefore, 
the dispersion of /3j n can disturb the Z-T> relation of a 
constant galactic age as shown in Figure 3. 

3. 3. Implication for Very Young Galaxies 

All our models considered in figure 6 predict that the 
dust-to-metal ratio in a galaxy younger than about 1 Gyr 
is very small, about one-fifth of the present-day Galactic 
value (dust/metal < 0.1). Because of the instantaneous 
recycling approximation, the accuracy of the model re- 
sult may not be so good for such very young galaxies. 
Nevertheless, we discuss what our model implies against 
very young galaxies whose age is less than 1 Gyr. 

First, we compare our model prediction with observa- 
tions of a nearby dwarf galaxy, SBS 0335—052. This dwarf 
galaxy has an age of less than about 100 Myr (Izotov, 
et al. 1997). In figures 8 and 9, we show such com- 
parisons. The allowed range of the metallicity is taken 
from Izotov, et al. (1997). That of the dust-to-gas ratio is 
based on the infrared spectral energy distribution (Dale, 
et al. 2001) and the hydrogen Brackett lines ratio (Hunt, 
et al. 2001). For theoretical predictions, we assume t$f = 1 
(thick curves) and 5 (thin curves) Gyr. As can be seen 
in figure 8, the effect of the metallicity in the infalling 
gas, which is assumed to be 0.01-Zq, appears in the solid 
(r/ = 0.5) and dashed (n = 0.3) curves. The dash-dotted 
and dotted curves are the cases of zero metal infall and no 
infall, respectively. From these figures, we find that our 
model reproduces the observed metallicity, dust-to-gas ra- 
tio, and very young age quite well if we assume tsf = 1 




Fig. 8. Metallicity evolution at a galactic age of less than 1 
Gyr. The solid curves are the standard case. The dashed, 
dotted, and dash-dotted curves are the cases of r\ = 0.3, no 
infall, and Z; n = 0, respectively. The assumed Tgp is 1 Gyr 
and 5 Gyr for the thick and thin curves, respectively. The 
shaded areas are the areas indicated from observations of a 
nearby dwarf galaxy, SBS 0335—052, and the Lyman break 
galaxies at z ~ 3, as denoted in the panel. The vertical axis 
is normalized by Zq. 

Gyr. Such a short time-scale of the star formation, i.e., a 
starburst like time-scale, is reasonable for SBS 0335—052 
because the galaxy is really starbursting now. 

Next, we move on a comparison with the Lyman break 
galaxies (LBGs) at z ~ 3. The LBGs seem to have an age 
from about 100 Myr to 1 Gyr (e.g., Shapley, et al. 2001) 
and a metallicity around one-third of the Solar value (e.g., 
Pettini, et al. 2001). The dust-to-gas ratio in the LBGs is 
still quite uncertain. We compare only the metallicity of 
the LBGs with our model, figure 8 shows that our models 
with Tgp = 1 Gyr reproduce the observed metallicity and 
age of the LBGs very well. In figure 9, we present the 
expected range of the dust-to-gas ratio in the LBGs (the 
thin shaded area); 0.01-0.1 of the dust-to-gas ratio of our 
Galaxy. Our model predicts that the dust-to-metal ratio 
in the LBGs is ~ 0.1 or less. This prediction should be 
tested in the future. 

For very young galaxies, our model predicts a lower 
dust-to-metal ratio than that in Edmunds (2001) and 
Morgan, Edmunds (2003). This is caused by the differ- 
ent assumptions between both models, for example, a low 
or high accretion efficiency (sticking coefficient), effective 
or no SNc destruction, with or without galactic gas infall. 
However, our smaller dust-to-metal ratio does not predict 
more transparent young galaxies straightforwardly. This 
is because the amount of attenuation depends on the dust 
optical depth, which is strongly affected by the geometry 
of the dust distribution. If the dust grains are confined 
in a small large attenuation is expected, even in 

a small dust amount. Of course, the dumpiness of the 
distribution is also important. 



Fig. 9. Dust-to-gas ratio evolution at a galactic age of less 
than 1 Gyr. The thick solid, dashed, dash-dotted, and dot- 
ted curves are the cases of the standard, no SNe destruction, 
high accretion growth, and no infall with rgp = 1 Gyr, re- 
spectively. The thin solid and dotted curves are the standard 
and no infall with rgp = 5 Gyr, respectively. The effect of 
metallicity in the infalling gas is negligible in this plot. The 
thick shaded area is the area indicated from the observations 
of SBS 0335—052. The thin shaded area is the expected area 
of the Lyman break galaxies at z ~ 3. The vertical axis is 
normalized by the dust-to-gas ratio of our Galaxy. 

4. Conclusion 

We construct a simple evolutionary model of the 
amount of metal and dust in order to discuss the evo- 
lution of the dust-to-metal ratio in galaxies. Then, the 
achieved conclusions are summarized as follows: 

(1) We propose that the grain growth time-scale by the 
metal accretion normalized by the gas metallicity (r acCi o) 
does not depend on the star- formation history (tsf) sig- 
nificantly. If this is correct, the ratio of the time-scale 
of dust destruction by SNe (tsn) to that of the accre- 
tion growth depends on the star-formation history because 
only tsn depends on tsf- On the contrary, a constant ra- 
tio rsN/Tacc,o has been assumed in previous studies, such 
as Dwek (1998). 

(2) The evolutionary track in the plane of the metallic- 
ity (Z) and the dust-to-gas ratio (T>) depend on the ratio 
tsn/tuccO, or equivalently on the star-formation history, 
tsf- 

(3) The observed linearity of the Z-V relation for 
nearby spiral galaxies is interpreted as the sequence of 
a constant galactic age in our framework. 

(4) The dust-to-metal ratio can evolve significantly dur- 
ing the evolutionary time-scale of galaxies: the ratio may 
be 50%, 30% and 20% of the local value at z ~ 0.5, 1, 
and > 2, respectively. Therefore, we will observe Z-T> re- 
lations different from the local one for high-z (or young) 
galaxies. 

(5) The uncertainties concerning the time-scales of dust 
destruction by SNe and dust growth by accretion are still 
large, and different parameter choices result in different 
evolutions of the dust-to-metal ratio. Fortunately, we 
can put a further constraint on these time-scales by using 
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the observations of the dust-to-metal ratio in high-z (or 
young) galaxies. Therefore, such observations are strongly 
encouraged. 

(6) The effect of infall on evolutionary tracks in the Z-V 
plane is negligible. However, since the speed of chemical 
enrichment depends on the infall rate, the dispersion of 
the infall rate can disturb the Z-T> relation of a constant 
galactic age. 

(7) Our model shows a very good agreement with the 
observations of a nearby extremely young (age < 100 
Myr) dwarf galaxy, SBS 0335-052. We suggests that 
the Lyman break galaxies at z ~ 3 have a low dust-to- 
metal ratio of 0.1, or less. Although our model predicts 
a relatively smaller amount of dust in high-z galaxies, it 
does not mean the transparent high-z galaxies, because 
the dust opacity strongly depends on the geometry of the 
dust distribution. 

The author thanks the referee, Dr. M. G. Edmunds, 
for his great efforts to understand this work, and also 
thank Ryuko Hirata and Hideyuki Kamaya for continu- 
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